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An analysis is given of thermoelastic noise (thermal noise due to thermoelastic dissipation) in finite 
sized test masses of laser interferometer gravitational-wave detectors. Finite-size effects increase 
the thermoelastic noise by a modest amount; for example, for the sapphire test masses tentatively 
planned for LIGO-II and plausible beam-spot radii, the increase is < 10 per cent. As a side issue, 
errors are pointed out in the currently used formulas for conventional, homogeneous thermal noise 
(noise associated with dissipation which is homogeneous and described by an imaginary part of the 
Young's modulus) in finite sized test masses. Correction of these errors increases the homogeneous 
thermal noise by < 5 per cent for LIGO-lI-type configurations. 



PACS numbers: 04.80.Nn, 05.40.-a 



I. INTRODUCTION AND SUMMARY 

Braginsky, Gorodetsky and Vyatchanin ]jj (henceforth 
BGV) have recently shown that noise associated with 
thermoelastic dissipation in test masses will affect signif- 
icantly the performance of LIGO-II gravitational-wave 
detectors. f\ 

The BGV computation of thermoelastic noise was 
based on an idealization in which each test mass has an 
arbitrarily large radius and length compared to the size 
of the light's beam spot on the mirrored test-mass face. 
In this limiting case, BGV showed that the spectral den- 
sity Sh(f) of the thermoelastic gravitational- wave noise 
scales as the inverse cube of the beam-spot radius r , 
Sh oc so it is desirable to make r Q large. However, 

when r a is no longer small compared to the test-mass 
size, the BGV analysis breaks down. 

The principal purpose of this paper is to explore, quan- 
titatively, the sign and magnitude of that breakdown. 
As we shall see, that breakdown (i.e., finite size of the 
test masses) increases the thermoelastic noise; but for 
expected beam-spot radii (r < 3/10 the test-mass ra- 
dius a), the increase is modest (< 10 per cent). 

A second purpose of this paper is to show how the 
BGV analysis of thermoelastic noise can be simplified 
considerably; and (adapting techniques due to Bondu, 
Hello and Vinet || - henceforth BHV), to show how to 
generalize the BGV analysis to finite sized test masses. 

A third purpose is to point out and correct errors in the 
BHV analysis of conventional homogeneous thermal noise 
(noise associated with a small, position- independent, 
imaginary part of the Young's modulus) in finite sized 



* "LIGO-II" refers to the second generation of detectors in 
the Laser Interferometer Gravitational Wave Observatory M 
and in its international partners. 



test masses. The corrections of the BHV formulae in- 
crease homogeneous thermal noise by < 5 per cent for 
beam-spot radii < 3/10 the test-mass radius a, and thus 
are primarily of conceptual importance, not practical im- 
portance. 

In Sec. 0, we outline our method of computing ther- 
moelastic noise, in Sec. [II we use our method to verify 
the BGV result for thermoelastic noise in the limit of 



arbitrarily large test masses, in Sec. IV we compute the 
thermoelastic noise in finite sized test masses and esti- 
mate the accuracy of our analysis, in Sec. ^ we correct 
the errors in the BHV computation of conventional, ho- 
mogeneous thermal noise, and in Sec. |v| we make some 
concluding remarks. 



II. METHOD OF CALCULATION 

Our analysis of thermoelastic noise is a simplification 
of one of the procedures developed by BGV: Appendix C 
of Ref . [0 . The foundation of the analysis is Levin's Q] 
"direct" method of computing thermal noise (of which 
thermoelastic noise is a special case): 

Levin begins by noting that the gravitational- wave de- 
tector's laser beam reads out a difference of generalized 
positions q(t) of the detector's four test masses, with each 
q given by an average, over the beam spot's Gaussian 
power profile, of the normal displacement 8z = u z of the 
test-mass face: 



- — Sz(r, <f))rd(f)dr 



Sz(r, (f>) rdipdr 



(1) 



Here (r, (f>) are circular polar coordinates centered on the 
beam-spot center (which we presume to be at the center 
of the circular test-mass face), r is the radius at which 
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the spot's power flux has dropped to 1/e of its central 

2/2 

value, and a is the test-mass radius. (The factor e~ a < r ° 
must be <C 1 in order to keep diffraction losses small, so 

2 I 2 

we shall approximate f — e~ a ' r ° by unity throughout this 
paper.) Levin then appeals to a very general formulation 
of the fluctuation-dissipation theorem, due to Callan and 
Welton H , to show that the test-mass thermal noise can 
be computed by the following thought experiment: 

We imagine applying a sinusoidally oscillating pres- 
sure, 



P = F n 



2 / 2 

— — cos(W) 



(2) 



to one face of the test mass. Here F a is a constant force 
amplitude, u> = 2nf is the angular frequency at which 
one wants to know the spectral density of thermal noise, 
and the pressure distribution (||) has precisely the same 
spatial profile as that of the generalized coordinate q, 
whose thermal noise S q (f) one wishes to compute. 

The oscillating pressure P feeds energy into the test 
mass, where it gets dissipated by thermoelastic heat flow. 
We compute the rate of this energy dissipation, Wdiss, av- 
eraged over the period 2it/uj of the pressure oscillations |] 
Then the fluctuation-dissipation theorem states that the 
spectral density of the noise S q (f) is given by 



S g (f) 



8k B TW dit 
F?uj 2 



(3) 



(Eq. (2) of Ref. (U); here fee is Boltzman's constant The 
interferometer's gravitational- wave signal h(t) is the dif- 
ference of the generalized positions q of the four test 
masses, divided by the interferometer arm length L. Cor- 
respondingly the contribution of the test-mass thermoe- 
lastic noise to the gravitational- wave noise is 1/L 2 times 
the sum of S q (f) over the four test masses (which might 
have different beam-spot sizes and thus different noises): 



(4) 



A=l 



The rate Wdiss of thermoelastic dissipation is given by 
the following standard expression (first term of Eq. (35.1) 
of Landau and Lifshitz Kil, cited henceforth as LL): 



/ TdS \ 



-(VST) 2 rd<j)drdz) . (5) 



^It is here that our analysis is simpler than that of BGV. 
Instead of computing Wdiss and using Eq. Q for the thermal 
noise, BGV compute the imaginary part Ss(x) of the test- 
mass susceptibility \ (which is much harder to compute than 
Wdiss) and then evalute S q in terms of 3(x) [their Eq. (14)]. 



Here the integral is over the entire test-mass interior us- 
ing cylindrical coordinates; T is the unperturbed temper- 
ature of the test-mass material and ST is the temperature 
perturbation produced by the oscillating pressure; dS/dt 
is the rate of increase of the test mass's entropy due to 
the flux of heat — kVST flowing down the temperature 
gradient VST, k is the material's coefficient of thermal 
conductivity, and (...) denotes an average over the pres- 
sure's oscillation period 1// = 2tt/lu. (For conceptual 
clarity we explicitly write the average (...) throughout 
this paper, even though in practice it gives just a simple 
factor (cos 2 uit) = 1/2.) 

To compute the thermal noise, then, we must calculate 
the oscillating temperature perturbation ST(r, cj), z,t) in- 
side the test mass, perform the integral (|5|) over the test- 
mass interior and the time average to obtain the dissipa- 
tion rate, then plug that rate into Eq. (pi) and then Eq. 
(©■ 

The computation of the oscillating temperature per- 
turbation is made fairly simple by two well-justified ap- 
proximations [jlj: 

First: The radius and length of the test mass are 
a H ^ 14 cm and the speeds of sound in the test-mass 
material are c s ~ 5 km/s, so the time required for sound 
to travel across the test mass is r soun( j ~ 30 /is, which 
is ~ 300 times shorter than the gravitational-wave (and 
pressure-oscillation) period r gw = 1// ~ 0.01 s. This 
Tsound <C T gw means that we can approximate the oscilla- 
tions of stress and strain in the test mass, induced by the 
oscillating pressure P, as quasistatic. It seems reasonable 
to expect this approximation to produce a fractional er- 



^quasistatic 



/sound 



/ 



300 (6) 



in our final answer for the thermoelastic noise S q (f). 
Here 



1 



sound 



rsound min(a, H) 



30, 000Hz 



(7) 



for the currently contemplated LIGO-II test masses: sap- 
phire with a ~ H ~ 14 cm. 

Second: The timescale for diffusive heat flow to alter 
the temperature distribution, tt ~ Cvpr 2 /n ~ 100 s, is 
~ 10 4 times longer than the pressure-oscillation period 
T gw (here CV — 7.9 x 10 6 erg g^ 1 K _1 is the specific heat 
per unit mass at constant volume, p ~ 4.0 g/cm 3 is the 
density, r Q ~ 4 cm is the spot size and k ~ 4.0 x 10 6 erg 
cm s _1 K _1 is the thermal conductivity, and our val- 
ues are for a sapphire test mass). This tt S> T gw means 
that, when computing the oscillating temperature dis- 
tribution, we can approximate the oscillations of stress, 
strain and temperature as adiabatic (negligible diffusive 
heat flow). The only place that heat flow must be consid- 
ered is in the volume integral (||) for the dissipation. The 
dominant contribution to that volume integral will come 
from a region with radius ~ r Q and thickness ~ r a near 
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the beam spot. The region of the integral in which the 
adiabatic approximation breaks down is predominantly a 
thin "boundary layer" near the beam spot with radius r 
and thickness of order the distance that substantial heat 
can flow in a time ~ l/^gw = l//> he., thickness of order 



CvPf 



/100Hz „ 
0.4mmW — - — for sapphire. (8) 



This region of adiabatic breakdown is a fraction ~ 
'"heat/^o °f the region that contributes substantially to 
the integral, so we expect a fractional error 



Hicat n m 
^adiabatic ~ ~ V.Ul 



(9) 



in S q (f) due to breakdown of the adiabatic approxima- 
tion. 

The quasistatic approximation permits us, at any mo- 
ment of time t, to compute the test mass's internal dis- 
placement field u, and most importantly its expansion 



9 = V • u , 



(10) 



from the equations of static stress balance (Eq. (7.4) of 
LL 0) 



(1 - 2a)V 2 u + V(V • u) = 



(11) 



(where a is the Poisson ratio), with the boundary con- 
dition that the normal pressure on the test-mass face be 
P(r, t) [Eq. (||)] and that all other non-tangential stresses 
vanish at the test-mass surface. Once has been com- 
puted, we can evaluate the temperature perturbation ST 
from the law of adiabatic temperature change (Eq. (6.5) 
ofLL§) 



ST 



-a t ET 



C v p(l - 2a) 



e 



(12) 



here a; is the linear thermal expansion coefficient, E is 
Young's modulus and Cy is the specific heat per unit 
mass at constant volume. [] This temperature perturba- 
tion can then be plugged into Eq. (|H|) to obtain the dis- 
sipation Wdiss as an integral over the gradient of the ex- 
pansion 



Em V 



(V6) 2 rd(j)drdz 



(13) 



This Wdiss can be inserted into Eq. (||) to obtain the 
thcrmoelastic noise. 



III. INFINITE TEST MASSES 

A. Dissipation and noise computed via BGV 
techniques 

We illustrate the above computational procedure by 
using it to verify the BGV Q result for thermoelastic 
noise in the case where each test mass is arbitrarily large 
compared to the spot size. 

Following BGV, we approximate the test mass as 
an infinite half space. Then the solution to the qua- 
sistatic stress-balance equation ([ll]) is given by a Green's- 
function expression [LL Eq. (8.18) with F x = F y = 0, 
F z = P(r, </))], integrated over the surface of the test 
mass. Taking the divergence of that expression [or, equiv- 
alently, taking the divergence of Eq. (39) of BGV], we 
obtain the following equation for the pressure-induced 
expansion: 

(1 + <t)(1-2<t)F 
6 = „ „ cos(wi) 



X Z 



ir 2 r 2 E 
dx'dy' ■ 



-(x' +y' )/r 



[(x - x') 2 + (y— y') 2 + z 2 } 3 / 2 



(14) 



where we have converted from polar coordinates to Carte- 
sian coordinates. Following a clever procedure implicit in 
the BGV analysis [in going from their Eq. (39) to (40)], 
we insert into the integral ( p^) an integral of the Dirac 
delta function written as 



8{x - x' - x")dx" 
1 



2tt 



e *A*-*'-x") dkxdx " (15) 



and a similar expression for J 5(y — y' — y")dy" , and we 
rewrite x — x' and y — y' in the denominator as x" and 
y" , thereby obtaining a new version of (|l4j) with integrals 
over k x ,k v ,x' ,y', x" , y". The integrals over x',y', x", y" 
are then readily carried out analytically (they are well- 
known Fourier transforms), to yield Eq. (40) of BGV|] 

(1 + <7)(1-2<7)F 

8 = cos mt 



2n 2 E 



+ 00 



-klrl/4 e -k ± z e i(k x x+k y y) dkxdky ; ( 16 ) 



where k± = \/k 2 + k 2 . 



It is straightforward to take the gradient of this ex- 
pression, square it (with one term an integral over k x , k y 



*LL use the volumetric thermal expansion coefficient a = 
3a; and the specific heat per unit volume C v = pCv- 



§Note that our notation differs slightly from that of BGV: 
Our x is their z, our z is their x, and they have factored out 
the costirf, which they write as e 1 ^'. 
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and the other over k' x ,k'), and integrate over x and 
y (from — oo to +00) and over z (from to 00); the 
result is J (VO) 2 dxdydz expressed as an integral over 
x,y,z,k x ,k y ,k' x ,k' y . The integrals can be done easily, 
first over z to get l/(k± +k' 1 ), then over x and y to get 
Dirac delta functions, then over the k's. The result, when 
inserted into Eq. (IHf), is 



(1 



2 „2„3 



(17) 



By then inserting this into Eq. (||), we obtain the BGV 
result for the thermoelastic noise [their Eq. (12)] 



S ™(/) = 



8(1 + a) 2 nafk B T 2 



^irC^p 2 ^ 



(18) 



Here the superscript ITM means for an "infinite test 
mass." 



The nonzero components of the gradient of this expansion 
are 



dO _ 2fi 
dr A + /i J 
80 _ 1[i 
dz X + fij 



a(fc)e~ fcz Ji(fcr)fc 3 dfc , (23a) 
a(k)e- kz J (kr)k 3 dk . (23b) 



By squaring the gradient, integrating over the interior of 
the test mass, and using the relations 



J n (kr)J n (k'r)rdr = 



d(k-k') 
k 



(24) 



(which follow from the Fourier-Bessel integral), and by 
replacing the Lame coefficients by the Poisson ratio and 
Young's modulus [Eqs. (pl|)], and inserting the resulting 
J (y<3) 2 rd(f>drdz into expression (|l3|), we obtain the same 
result ( |Tt| ) as we got using BGV techniques. By inserting 
this into Eq. (13), we obtain the thermoelastic noise (p"8]). 



B. Derivation via BHV techniques 

Equation (0) for Wdi SS can also be derived in cylin- 
drical coordinates (r, z, (f>) using the techniques of BHV 
The displacement u has components [BHV Eqs. (5) 
and (6) with the denominator in (5) corrected from /1 to 
fx + A and with /? = a; see passage following BHV Eq. 
(8)] 

u r = f a(k) (l - ^±1^1 + kz \ e ~ kz ,h{kr)kdk , 



ti- = I a(k)(l + —^ h kz ) e~ kz J (kr)kdk 

lo \ X + p 

M0 = , 

where 



(19) 



a(k) 



A-Kjik 



- F„ cos ojt 



(20) 



[BHV Eq. (11), with the overall sign corrected from — 
to +, with w Q = \[2r cf. BHV Eq. (2), and with coswt 
inserted because our method of applying the fluctuation- 
dissipation theorem is dynamical while BHV's method is 
static]. In Eqs. ( |l9| ) the J„ are Bessel functions and A 
and n are the Lame coefficients (and \i is also the shear 
modulus), which are related to the Young's modulus E 
and the Poisson ratio a by 



A 



Ea 



E 



(1-2<t)(1 + ct) n 2(1 + a) 
The divergence of the displacement (|l^) is 
2/i 



(21) 



e = - 



^ + M Jo 



a(k)e~ kz Jo(kr)k 2 dk . (22) 



IV. FINITE SIZED TEST MASSES 
A. BHV Solution for Displacement 

Consider a finite sized, cylindrical test mass with ra- 
dius a and thickness H, and with the Gaussian shaped 
light spot centered on the cylinder's circular face. For 
this case, Bondu, Hello and Vinet (BHV) || have con- 
structed a rather accurate but approximate solution of 
the static elasticity equations. Unfortunately, their so- 
lution satisfies the wrong boundary conditions and thus 
must be corrected: 

The error arises when BHV expand the Gaussian- 
shaped pressure (g) as a sum over Bessel functions. BHV 
incorrectly omit a uniform-pressure term from the sum. 
As a result, the pressure that they imagine applying to 
the test-mass face [their Eq. (18)], 



-PbhvM = F a cos ^ PmJo{k m r) 



(25) 



[where Jo is the Bessel function of order zero, k m is re- 
lated to the m'th zero Q m of the order-one Bessel function 
J\{x) by k m — (, m /a, and p m are constant coefficients 
given below] , has a vanishing surface integral 



PBHV^nrdr = 



(26) 



In other words, their applied pressure (^5|) is equal to 
the desired pressure P(r) [Eq. @] minus an equal and 
opposite net force F Q cos(wi) applied uniformly over the 
test-mass face: 



Ft 



BHV 



(r) = P(r) -poF cos u>t ; 



(27) 



4 



Pa 



(28) 



2 / 2 

[Recall that we are approximating 1 — e~ a < r ° by unity; 
see discussion following Eq. (|l|).] 

It is evident, then, that to get the correct distribution 
of elastic displacement u inside the test mass, we must 
add to the BHV displacement a correction. This correc- 
tion is the displacement caused by the spatially uniform 
pressure poF a cos cot on the test-mass face. That uniform 
pressure causes the test mass to accelerate with acceler- 
ation a = [(F cos wt) / M]e z , where M — ira 2 Hp is the 
mass of the test mass and p is its density. In the reference 
frame of the accelerating test mass, all parts of the test 
mass feel a "gravitational" acceleration ge z equal and 
opposite to a, i.e. g = — (F a cosut)/M, (which can be 
treated as quasistatic, though it oscillates at frequency 
lj). Thus, the displacement is the same as would occur 
if the test mass were to reside in the gravitational field 
ge z with a uniform pressure on its face counteracting the 
force of gravity. The solution for this displacement is 
given by LL [|| (Problem 1, page 18).** Translating into 
our notation and converting from the Young's modulus 
and Poisson ratio to the Lame coefficients via Eq. (21), 
we obtain: 



Su r 



Xp r 



F a cos tot 2/i(3A + 2/i) 
8u z \p r 2 



z 



(29a) 



= (A + n)p a 

F cos ut ~ ApH(3X + 2/z) /i(3A + 2p) 



2H 



(29b) 



The total corrected displacement, in cylindrical coor- 
dinates, is 

u r = HV + 5u r , u z = v% HV + 8u z , U(P = , (30) 

where u? HV is the BHV displacement [their Eqs. (15) 
plus (25) and (17) plus (26)]: 



«« HV (r,z) A + 2/i 



F cos ujt 2/i(3A + 2/i) 

OO 

+ ^ An(z)Ji(k m r) , 

m— 1 

^ HV (r,z) _ A 



(c r + cirz) 



(31a) 



F a cos tot /i(3A + 2/i) 



c z - 



**LL seek to solve a problem in which (in the presence of the 
uniform gravitational acceleration), instead of having a uni- 
form pressure applied to the face of the cylindrical test mass, 
the face has vanishing displacement. Their solution actually 
satisfies our desired boundary conditions but not theirs; there- 
fore, they comment on it being inaccurate near the test-mass 
face. For our problem it is accurate. 



A + 2/i 
4^(3A + 2/i) 



c\r 2 + ^ B m (z)J (k ri 



,BHV 



(31b) 
(31c) 

Here the coefficients Co and c\ are [equations following 
Eqs. (24) and (26) of BHV] 



(r,z) = 



a 2 ^ J (( m )p r , 



rn — 1 



Sti 



Cl 



-2c 
H 



(32) 



and A m and B m are the following functions of z [Eqs 
(19) and (20) of BHV] 

~t~ ^m^ 

k m z A + p i _f. 



B m {z) 



2 X + 2p 
A + 3/i 



2(A + 2/x) 
X + 3p 



2(A + 2/i) 
k rn z A -+- p 



Pm - 5, 



-k m z 



(33) 



2 A + 2/i f """ > ' (34) 

where a m , /3 m , 7 m and S m are constants given by [Eqs. 
(21)-(24) of BHV]: 

Q m = exp(-2fc m iJ) 

Pm 

(A + 2/i) 1 - Q m + 2k m HQ m 



7n 



k m p(X + p) (1 - Q m ) 2 - Ak 2 m H 2 Q m 
p m (X + 2p)Q m 1 - Q m + 2k m H 
k m p(X +n) (1 - Q m ) 2 - ±k 2 n H 2 Q m 

Pm 

2k m p(X + p) 
[2k 2 n H 2 (X + p) + 2pk m H]Q m + /i(l - Q m ) 
(1 - Qm) 2 - 4k 2 m H 2 Q m 



(35a) 
(35b) 



2k m p(X + p) 

2k 2 m H 2 (X + /i) - 2pk m H - /i(l - Q m ) 
X (1 - Qm) 2 ~ 4k 2 m H 2 Q m 

with [equation following Eq. (18) in BHV] 

2 f a e~ r2 / r ° 
Pm = 7 7277 \ / 5— Jo(k m r) r dr . 



(35c) 



(35d) 



(35e) 



a 2 Jfi(( m ) Jo nr 2 



(36) 



2/2 

In the spirit of our approximating 1 — e~ a < r ° by unity 
[discussion following Eq. ([[])], BHV suggest approximat- 
ing the upper limit of this integral by 00; the integral can 
then be done analytically, giving [equation preceding Eq. 
(19) of BHV] 



Pm = 



exp (-fc^rg/4) 

1,2 J Q (Cm) 



(37) 



This is a good approximation to the exact formula ( |3q ) 
for small m (which turn out to give the dominant con- 
tribution to the noise), but for large m it can severely 
underestimate p m . 
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B. Expansion and the integral of its squared gradient 

It is straightforward to compute the expansion = 
V • u and the components of its gradient from expressions 
(29), (fl|) and @); the results are 

Q(r,z) _ po ( z\ x 2(c + ciz) 



F cosujt 3A + 2 ( u\ HI 3A + 2/i 

+ Y / [k m A m {z) + B' m {z)] J (k m r) , (38) 

m 

and 

® Q/dr = - y] k m [k m A m (z)+B' m (z)] Ji [k m r), (39a) 

i- q COS LUL 

m 

<MI ' h ' + B'L{z)] Mk m r) , 

tn 

(39b) 



F cos wt 3\ + 2[i 



where the primes denote derivatives with respect to z and 
the coefficient c\ is 



Ci = Ci 



Po 
2H 



(40) 



Using the (nonstandard) orthogonality relations 



r J\{k m r) Ji(k n r)dr = — Jo(Cm) $mn , (41) 
a 2 

r J (k m r) J Q (k n r) dr = — J (( m )5 mn , (42) 



r J (k m r) dr = 



(43) 



the volume integral of (VO) 2 can be evaluated analyti- 
cally. The result, after some algebra and after averaging 
cos 2 Lot to 1/2, is 

1 / [rim , , \ 2ira 2 £ 2 H 



E 



^mPm(l Qm)Jo(Cm) 



(44) 



2 < A + ' 1 ) 2 „fi[(i-W 2 -^TO 

X [(1 - g m ) 2 (l + Qm) + 8Hk m Qm(l - Qra) 

+4H 2 kf n Q m (l + Q m )] . 



C. Thermoelastic noise 

Inserting Eq. into Eq. ( [To] ) and then into Eq. (||), 
and using Eqs. ( pi|) for the Lame coefficients, we obtain 
for the spectral density of thermoelastic noise in a finite 
sized test mass: 



rtFTM _ n 2 oITM 



(45) 



Here S^™ is the BGV result @ for the spectral density 
for an infinite test mass, and C FTM is the following finite- 
test-mass correction to the spectral density: 



r 2 

°FTM 



E 



a 5 H£ 2 
(X k m p m (l — Qro)^o(Cm) 

i [(1 ~ Qm) 2 ~ AH 2 k 2 m Q m f 

[(l-Qm) 2 (l+Qm)+8Hk„ 

+4H 2 k 2 m Q m {l + Q m )] } . 



(46) 



The square root, CftMj of this finite-test-mass correc- 
tion is plotted in Fig. |l| as a function of the test-mass 
thickness H and radius a measured in units of the beam- 



spot radius r Q . [One can easily show from Eq. (46) that 
Cftm depends on H , a and r Q only through the dimen- 
sionless ratios H/r and a/r a , as must be the case on 
dimensional grounds.] Notice that the noise is larger, 
at fixed r a , for large-a, small- H test masses (thin disks) 
than for small-a, large- H test masses (long cylinders). 
However, for plausible parameters the difference is only 
a few tens of per cent. The reason for the greater noise 
in a thin disk is that it experiences greater deformation, 
when a force acts at the center of its face, than does a 
long cylinder, and thus the integral (jl^), to which the 
noise is proportional, is larger. (See, e.g., Sec. 12 of ||, 
or Sec. 305 of §.) 
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FIG. 1. Contour plot of th e corr ection Cftm to the ther- 
moelastic amplitude noise \/ S„(f) due to the finite size of 
the test mass [Eqs. (j^|) and (^6|)]. This correction is shown 
as a function of the test-mass radius a and thickness H, both 
measured in units of the beam-spot radius r (the radius at 
which the light beam's energy flux has dropped to 1/e of its 
central value.) 

The current "straw- man" ("reference") design for 
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LIGO-II includes sapphire test masses with a = 14 cm 
and H = 12.2 cm. In Fig. ||we plot the finite-test-mass 
correction Cftm as a function of beam-spot radius r a 
(in centimeters) for such test masses (for which we use 
the BGV values of the parameters a = 5.0 x 10~ 6 K _1 , 
k = 4.0 x 10 6 erg K _1 cm -1 s -1 , p = 4.0 g/cm 3 , 
C v = 7.9 x 10 6 erg g" 1 K -1 , E = 4 x 10 12 erg/cm 3 , 
(7 = 0.29). Although we continue our plot up to r a = 6 
cm, it may be impractical or undesirable to operate with 
r much larger than 4 cm. Two reasons for this are: (i) 
Each time the light beam encounters a test mass, a frac- 

2/2 

tion ~ e~ a > r ° of its power is lost around the test-mass 
sides ("diffraction losses"); keeping this below ~ 10 ppm 
requires r a < 4 cm. (ii) There are practical limitations 
R < 50 km on the radii of curvature of the test-mass 
mirrors; if the beam waist is half way between the mir- 
rors of an arm's optical cavity so the spot sizes r D are 
the same on the two mirrors, and if R is significantly 
larger than the arm length L — 4 km, then the spot sizes 
are r ~ (X 2 LR/8n 2 ) 1 ^ 4 (where A = 1.06/im is the light 
wavelength), so R < 50 km requires r < 4 cm. 

For the plausible range r a < 4 cm, Fig. || shows that 
the finite-test-mass correction to the amplitude noise is 
< 10 per cent. 



fractional error e a diaba 



0.01 in S q [Eq. ©], and the 




3 

r ,cm 

FIG. 2. Correction Cftm to the thermoelastic amplitude 
noise as a function of the beam-spot radius r , for test 
masses with the parameters currently being contemplated for 
LIGO-II: 30 kg sapphire cylinders with radius a = 14 cm and 
thickness H = 12.2 cm. 



D. Errors in our analysis 

There are three significant sources of error in our anal- 
ysis. We expect them to produce a net error in Cftm 

and thence in the test-mass noise J S q ™ that is < 1 per 

cent, for the expected LIGO-II parameter regime (a ~ 14 
cm, H ~ 12 cm, r a < 4 cm). More specifically: 

One error source is the quasistatic approximation. We 
have already estimated this as producing a fractional er- 
ror £q U asistatic < 0.003 in S h [Eq. @], and the error in 
^fS~ q will be half this, < 0.0015. 

The second error source is the adiabatic approxima- 
tion. We have already estimated that this produces a 



error m 



will be half this, < 0.005. 



The third error source is one that we have not dis- 



cussed: A failure of the elastic displacement (31) to sat 
isfy the boundary condition T, 
cylindrical sides, r = a 



on the test mass's 
As was discussed by BHV [|, 
the Co and ci terms in the displacement (31 ) are a correc- 
tion to the leading-order displacement, designed to im- 
prove the satisfaction of the T rr (a) = boundary condi- 
tion. We shall refer to these terms as the "Saint- Venant 
correction" In our final answer for S q (f) [Eqs. ( |45| ) 
and (|iq)], this Saint- Venant correction makes a fractional 
contribution < 6 per cent, for LIGO-II type test masses 
and plausible beam radii r Q < 4 cm. The rms value 
of T rr {a) with the Saint- Venant correction included is 
smaller than that without the Saint- Venant correction 
by about a factor 3, so it is reasonable to expect that the 
remaining error in S q (f) due to T rr (a) =/= is < 1/3 of 
the Saint- Venant correction, i.e., a remaining fractional 
error 



1 



esv < 77 x 0.06 = 0.02 



(47) 



% will be half this, < 0.01 



The fractional error in 
which is larger than the other two errors 

Combining these three errors in quadrature, we expect 

our formulas for 

of magnitude 



S q ™ to make a net fractional error 



Ar SJS™ 

(ICFTM V 1 



Cftm /^ftm 



< 0.01 



(48) 



for LIGO-II-type test masses and beam-spot radii r a < 4 
cm. 



V. CONVENTIONAL THERMAL NOISE 

Because of the boundary-condition error that BHV 
make in solving the elasticity equations (and because of 
an additional algebraic error discussed below), their re- 
sult for the conventional thermal noise must be corrected. 

The conventional thermal noise is given by Levin's for- 
mula (^) with Wdiss the time-averaged dissipation pro- 
duced by an imaginary part ^s{E) — &(lu)E of the 
Young's modulus: 

W diss = 2u/$(w)(lT) 

= w$(w) J (AO 2 + 2nSijSij)rdrd(j)dz . (49) 

Here (U) is the time-averaged elastic energy, SijSij is the 
square of the strain associated with the displacement u, 
there is an implied sum over i and j, and the integral is 
over the test-mass interior; cf. Eq. (12) of Ref. B. 
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The expansion O is given by Eq. (pq), and the com- 
ponents of the strain on the spherical, orthonormal basis 
eV, e^, e z are readily computable from the displacement 
©, (29), fU) via [Eqs. (A.1)-(A.4) of BHV] 



r 



(50) 



where commas denote partial derivatives. By evaluating 
these strain components, inserting them and the expan- 
sion ( |38| ) into Eq. ([49]), averaging over time, integrating 
over the test mass, and reexpressing the Lame cocfficnts 
in terms of the Young's modulus and Poisson ratio, we 
obtain 



W diss = cj$(u;) (U„ + AU) . 
Here U is given by 

(1-^3 » p^JliU) 



u n = 



E^ 



E f— ' Cm 

m— 1 

with [equation following Eq. (29) of BHV] 
1 - Qi + ik m HQ m 



while AZ7 is 



AC/ = 



(1 - Q m ) 2 - 4k 2 m H 2 Q r . 



(51) 



(52) 



(53) 



6ttH 3 E 



[ir 2 H 4 pl + 12TrH 2 <7p s + 72(1 - o)s 2 } , 

(54) 



with 



! E 

m— 1 



PmJo(Cm) 

c 2 



(55) 



When the approximation ( |37| ) is made for p m , [/ G takes 
the form given by BHV [their Eq. (30)] 



1 z 

m— 1 

and 5 takes the form 



exp(-C r 2 »rg/2 fl 2 ) 

Cm^o(Cm) 2 



E 



exp(-Cr 2 /4a 2 ) 

C m ^0(Cm) 



(57) 



The approximations ( pq ) and ( j57j ) are rather good for 
realistic parameter values, despite the fact that for large 
m Eq. (p7| ) is a very poor approximation to p m , because 
large m make small contributions to U and s. 

Equations Q) and ( p7| ) for AJ7 differ from Eq. (31) of 
BHV for two reasons: (!) BHV used the wrong bound- 
ary c onditions at the test-mass face [see beginning of Sec. 
[V A above] ; correcting this leads to all the terms in Eq. 



(p4j) involving p a . (ii) BHV seem to have made an alge- 
braic error: Eqs. ( p4\ ) and ( |57| ) should agree with BHV 
Eq. (31) when p is set to zero, but they do not; it might 
be that BHV accidentally omitted the term or the 
S 2 r term when evaluating Eq. (|49|). 

Inserting Eq. ( |49| ) into Eq. (||), we obtain the BHV ex- 
pression for the conventional thermal noise [their equa- 
tion following Eq. (31)] 



cFTM (/) = 8fcBT $M([/o + Ac0 



(58) 



and 



where (to reiterate) U is given by Eqs. ( p2 f ) [o r ( p6[ ) 
(||), while AU is given by Eqs. (fj) and® [or @]. 

If the test mass is infinite in size, then the conventional 
thermal noise has the following form, derived by BHV 
[their Eq. (14) with w = \^2r a , which differs from the 
formula derived earlier by Levin B — his Eq. (2)]: 



S: 



ITM 



4fc R T 1 - a 2 



2nEr 



$(w) . 



(59) 



As for thermoelastic noise, we define a finite-test-mass 
2 the ra 
to that (t 



correction Cp TM to be the ratio of the finite-test-mass 



spectral density (| 
mass: 



for the infinite test 



Q 

r<i _ i 

U FTM 



FTM 



^ITM 



(60) 



We plot the square root of this correction (i.e., the 
amplitude-noise correction) as a function of beam-spot 
radius r in Fig. |3| for a LIGO-II type test mass. We 
show there two curves, Cftm as given by the BHV for- 
mulae, and as given by our corrected formulae. Note that 
the BHV errors have only a small influence: their noise 
was too low by a factor < 5 per cent when r a < 4 cm. 
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FIG. 3. Cftm = V ' S «™ I V ' S l™ > the rmite-test-mass 
correction to the conventional, homogenous thermal noise, as 
a function of the beam-spot radius r , for test masses with 
the parameters currently being contemplated for LIGO-II: 30 
kg sapphire cylinders with radius a = 14 cm and thickness 
H — 12.2 cm. The curve labeled "BHV" is the result derived 
in Ref. || [their Eqs. (29), (28), and equation following (31)]; 
the curve labeled "corrected" is our corrected version of their 
result [our Eqs. @, @ and 
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VI. CONCLUSION 

In this paper we have sketched a fairly simple method 
of analyzing thermoelastic thermal noise in interferomet- 
ric detectors, we have used that method to derive for- 
mulas for the noise in cylindrical test masses with finite 
radius, thickness, and beam spots, and we have corrected 
the corresponding finite test-mass formulas for conven- 
tional thermal noise. Our formulas should be useful in 
optimizing the test-mass designs for interferometric grav- 
itational wave detectors. 

Because thermoelastic noise arises from physical pro- 
cesses associated with ordinary thermal fluctuations, 
thermal conductivity and thermal expansion, and is not 
influenced by "dirty" processes such as lattice defects and 
impurities (except through the easily measured conduc- 
tivity and expansion), the predictions for thermoelastic 
noise should be very reliable. Nevertheless, experimental 
tests of the theory would be useful and are being planned. 

Other forms of thermal noise do rely in crucial, ill- 
understood ways on dirty processes and thus are far less 
reliably understood than thermoelastic noise. This is es- 
pecially the case of thermal noise associated with (in- 
homogeneous) dissipation in and just beneath the test 
mass's dielectric-mirror coatings [for which Levin Gj pre- 
dicts, in the infinite-test-mass limit, a dependence S q cx 
X/Vg on beam-spot radius, compared to S q cx 1/r™ for 
thermoelastic noise and S q cx l/r for conventional, ho- 
mogeneous thermal noise] . Detailed experimental studies 
of these other forms of thermal noise are much needed as 
part of the R&D for interferometric gravitational-wave 
detectors, and are being planned. 

In some of the planned experiments, very small beam 
radii r Q and/or high frequencies / may be used. For 



„ < heat 
' o ^ I o 



CvPf 



/ 100Hz . , 

0.4mmW — - — for sapphire, (61) 



the adiabatic approximation breaks down seriously [cf. 
Eq. ([)]) and associated discussion] and our analysis of 
thermoelastic noise must be redone taking account of the 
diffusive redistribution of temperature during the elastic 
oscillations. Some foundations for doing this have been 
laid by BGV E|. For frequencies 



/>/, 



sound 



min(a, H) 
~ 10 4 Hz 



10cm 



min(a, H) 



for sapphire (62) 



(where c s is the sound speed), the quasistatic approxi- 
mation breaks down seriously [cf. Eq. (^]) and associated 
discussion], and our analysis must be redone taking ac- 
count of the finite propagation speed of the test mass's 
elastic deformations. 



After completing our analysis of thermoelastic noise 
in finite sized test masses, we learned that Sergey Vy- 
atchanin Q has been carrying out an analysis of this 
same issue, but using somewhat different techniques. In 
writing the final version of this paper, we have benefitted 
from email exchanges with him. 
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